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While the mechanical distortions change the electronic properties of graphene significantly, the 
effects of electronic manipulation on its mechanical properties have not been known. Using first- 
principles calculation methods, we show that, when graphene expands isotropically under equibiaxial 
strain, both the electron and hole doping can maintain or improve its ideal strength slightly and 
enhance the critical breaking strain dramatically. Contrary to the isotropic expansions, the electron 
doping decreases the ideal strength as well as critical strain of uniaxially strained graphene while 
the hole doping increases the both. Distinct failure mechanisms depending on type of strains are 
shown to be origins of the different doping induced mechanical stabilities. Our findings may resolve 
a contradiction between recent experimental and theoretical results on the strength of graphene. 

PACS numbers: 62.20.-x, 63.20.kd, 63.22. Rc, 81.40.Jj 



A defect-free infinite crystal becomes mechanically un- 
stable at a stress with a corresponding critical strain. 
Such a stress is the ideal strength of a material which 
is an inherent property of a given atomic and elec- 
tronic structure as well as a natural upper bound on its 
strength [U Hj]. While the variation of electron density 
can alter ideal strength of some materials 0], those ma- 
nipulations usually involve stoichiometric change which 
may alter the nature of chemical bonding between atoms 
in crystal. In this regard, graphene, a two-dimensional 
crystal with atomic thickness, is an ideal material to in- 
vestigate electronic control of ideal strength because it 
allows very high level of charge doping without any sacri- 
fice of atomic integrity through either field-effect transis- 
tor setups [4, 5] or alkali metal depositions [(| . Although 
there have been many studies about effects of mechani- 
cal perturbation on electron physics in graphene fj- 12|. 
the effect of electron density variation on the mechanical 
properties of graphene has not been known yet. 

Graphene is also known to be the strongest mate- 
rial [13| and can maintain its pristine strength regard- 
less of presence of large grain boundaries [l4|. A re- 
cent nanoindentation experiment [131 ] reports that the 
intrinsic strength of graphene is 42 N m _1 at the nom- 
inal equibiaxial breaking strain of e max = 0.225. This 
was followed by a theoretical calculation [ll| showing 
£max = 0.15 which is quite smaller than the measured 
value. It is unusual because usual materials under strain 
typically fail before they reach the ideal strength due to 
various reasons. Therefore considerations on the other 
external effects such as a doping are called for to un- 
derstand the ideal strength of graphene further since 
graphene samples in experiments are typically doped 
chemically or electronically. 

In this paper, we report a theoretical study show- 
ing that graphene under equibiaxial strain becomes to 
be stronger as doping increases regardless of the charge 
character of doping, electron or hole. Though the ideal 
strength under equibiaxial strain increases slightly by 



1 ~ 6 % depending on doping type and amounts, the 
corresponding critical breaking strain increases dramat- 
ically by ~ 20 % within the doping level of ~ 3 x 
10 14 /cm 2 which is accessible in experiments [1,0]. We 
show that the failure of isotropically expanded graphene 
through the A[ phonon mode softening [l5| or the two- 
dimensional Peierls instability [l6j] is overcome by dop- 
ing because the Fermi level of doped graphene lies out- 
side the energy gap associated with the A[ phonon. On 
the contrary, under uniaxial strain invoking elastic fail- 
ure of graphene [Hj], we find that electron doping weak- 
ens graphene by 5 ~ 7% while hole doping strengthens 
one by 3 ~ 6%. The corresponding breaking strains also 
change a lot by 37% as doping swept from 0.05 holes to 
0.05 electrons per carbon atom. Because the elastic fail- 
ure under uniaxial strain is associated with occupation 
of the er* band [13], the asymmetric dependence of ideal 
strength on electron and hole doping can be explained in 
term of doping induced occupation (electron doping) or 
depopulaton (hole doping) of the a* band of uniaxially 
strained graphene. 



Our first-principles calculations 17j are carried out 
with plane wave basis and norm-conserving pseudopo- 
tentials [HI- The local density approximation (LDA) is 
used for the exchange-correlation functional [lj| ■ Phonon 
spectrum is also calculated using density functional per- 
turbation theory 17 , 2(| • The energy cutoff for expansion 
of the plane wave basis set is 80 Ry. The /c-point grid of 
48 x 48 x 1 is used for the electronic structure and ionic 
structure relaxation calculations. Adjacent graphene lay- 
ers in neighboring supercells are 8 A away to avoid spuri- 
ous interactions between them. Hole and electron doping 
is simulated within the rigid band approximation which 
is successfully used in recent studies on Kohn anomaly 
of doped carbon nanotubes 21| and doped graphane [22j 
respectively. For phonon dispersion, 12 x 12 x 1 g-point 
sampling is used. Near the Kohn anomaly where the 
change of phonon frequency is very abrupt in a narrow 
momentum window, the explicit phonon calculations are 
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FIG. 1: (Color online) (a) Atomic displacements (red ar- 
rows) associated with the A[ phonon mode with an amplitude 
S a- The thick blue vectors indicate the Tf-point supercell of 
graphene. (b) The total energy change as a function of 8a un- 
der various equibiaxial strains (e). From top to bottom lines, 
e = 0.12, 0.14, 0.15, and 0.16. (b) and (c) show the same 
curves as in (a) at a fixed strain, e = 0.16, with changing 
electron and hole doping from top to bottom, n e (rih) = 0.04, 
0.02, 0.01, 0.00 respectively. 
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FIG. 2: (Color online) (a) and (b) show the stress-strain re- 
lationship with different electron and hole doping levels re- 
spectively. The red crosses are for undoped graphene for 
both (a) and (b). The equibiaxial strain is defined by <r c b = 
\J a% + <jy/y/2. In (a), n e = 0.013 (open green circles) and 
0.04 (solid blue). In (b), n h = 0.01 (open green) and 0.04 
(solid blue). The change of (c) ideal strength (<5o"id ea i) and 
(d) corresponding critical breaking strain (<5e C ritical) of doped 
graphene with respect to those of undoped one in percent- 
age as a function of n e (n^). In both plots, the empty (filled) 
rectangles are for electron (hole) doped graphene. 



additionally done along the high symmetry paths to avoid 
the inherent errors in interpolated phonon dispersions. 

When graphene expands isotropically, the A[ phonon 
mode (hexagonal breathing mode) at the K point soft- 
ens and then becomes to be unstable at a certain critical 
equibiaxial strain ■ The structures in Fig. Q] (a) show 
schematically how the crystal structure of graphene is de- 
formed when the A[ phonon is excited. Here, the strain 
is defined by e = (a — ao)/arj, where a and ao are the 
lattice constant of equibiaxially strained graphene and 
that of pristine graphene respectively. In Fig. [T](b), the 
total energy change as a function of the structural defor- 
mation of A[ phonon mode shows that, for e > 0.15, the 
hexagonally deformed structures are energetically more 
stable than the symmetric one 15]. For these cacula- 
tions we use the -fT-point supercell in which six carbon 
atoms are included as shown in Fig. 1(a) [15j . This can 
be understood as two-dimensional extension of Peierls 
instability 16]. A half- filled pristine graphene has its 



Fermi level at the Dirac point and the excitaion of A[ 
phonon mo de gene rates a dynamical energy gap at the 
Dirac point |16l. [23 1 realizing the Kohn anomaly at the K- 
point [23l . |24| . As graphene expands over critical equibi- 



axial strain, the energy gain through formation of static 
energy gap at the JsT-point overcomes the energy cost 
from ion-ion interaction so that the hexagonal network 
of graphene breakdowns 3, T^. Therefore, if the Fermi 
level shifts away from the Dirac point by doping, the 
Peierls instability will be avoided by hardening the A[ 
phonon mode so that graphene becomes more resistive 
against the applied equibiaxial strain. 

To confirm the doping related stability, we calculate 
the total energy of doped graphene with the A[ phonon 
mode. As discussed, Fig.[T](b) shows the lattice instabil- 
ity of undoped graphene over a critical equibiaxial strain 
(e > 0.15). We have done the same calculations for doped 
graphene at various doping levels with a fixed overcriti- 
cal strain of e = 0.16 [Figs. [T] (c) and (d)]. The electron 
(hole) doping level, n e (n^), is defined by the number of 
additional electrons (holes) per a carbon atom. We note 
that n e = 0.01 corresponds to 3.848 x 10 13 electron/cm 2 . 
As shown in Figs. [T] (c) and (d), the electron and hole 
dopings recover the stability of A[ phonon mode even at 
the strain of e = 0.16 under which pristine graphene is 
already broken. 

To obtain complete ideal strengths and corresponding 
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FIG. 3: (Color online) Phonon dispersion of (a) undoped and (b) electron-doped graphene (n e = 0.04). In both cases, red 
(blue) curves are for graphene with (without) equibiaxial strain ((a) e = 0.148 and (b) e = 0.167). The black arrows indicate 
the Kohn anomaly. The inset for each plot shows the vectors (black arrows) connecting Fermi surfaces (red solid lines) inducing 
the Kohn anomaly. Here we do not show irrelevant out-of-plane phonon mode spectrum, (c) <5eridcai (empty rectangles) in the 
left ordinate and <5e C riticai (filled rectangles) in the right one obtained from phonon dispersions as function of n e (abscissa) . 



critical breaking strains as function of doping, the stress- 
strain calculations are needed because the total energy 
calculations only indicates structural instability. In the 
stress-strain relationship, the maximum stress at the crit- 
ical breaking strain is the ideal strength of doped and un- 
doped graphene. Figs. [2] (a) and (b) show representative 
stress-strain curves of strained graphene with the electron 
and hole doping respectively. It is noted that the stress 
maximum corresponding to the ideal strength shifts to- 
ward larger strain as the doping level increases indepen- 
dent of the doping type. We obtain the ideal strength 
and the corresponding critical breaking strain as a func- 
tion of doping by taking the stress maximum values in 
Figs. [2] (a) and (b). The ideal strength (cidcai) of undoped 
graphene is found to be 44.6 N/m agre eing well with a 
recent experiment [l3| and calculation [l5| The critical 
breaking strain, however, is e = 0.148 much smaller than 
the measured value 15j, Il6 |. As n e increases from 0.00 
to 0.04, Cidoai increases by ~ 2% (n e — 0.02) and then 
decreases. As rih increases from 0.00 to 0.04, ordeal in- 
creases monotonically by ~ 6% [Figure 2(c)]. Contrary 
to that, the critical breaking strain (e cr jti C al) shows a dra- 
matic increase by ~ 18% (e cr iti ca i ~ 0.18) as both n e and 
nh increase to 0.04 [Figure 2(d)]. 

We have further calculated phonon dispersions in or- 
der to see how the phonon softening at the A"-point is af- 
fected by the doping. Softening of phonon to the negative 
phonon frequency may precede or indicate the instability 
of crystal structure [25(. Fig. [3] (a) shows the phonon 
dispersion of graphene with e = (red) and e = 0.148 
(blue) agreeing well with a previous study [13]. The fre- 
quency dip at the if-point associated with the A[ mode 
is the Kohn anomaly 2^, 24 1. We note that the LDA is 
not quantitatively perfect to describe the anomaly [26| 
but is sufficient to describe the softening within a cer- 
tain error bar 15] ■ The Kohn anomaly occurs only for q 
when two electronic states at ki and ki + q are on the 



Fermi surface 2J, |27| . In undoped graphene, the Fermi 
surface is at the Dirac points located at K and K' = 2K 
so that q = (r-point) and q = K (if-point) are al- 
lowed for the anomaly. The inset in Fig. [3] (a) shows the 
corresponding vector K in the BZ. As graphene doped, 
the momentum at which the Kohn anomaly occurs shifts 
from the if -point and splitts into two points. The A\ 
phonon hardens also because of change of the Fermi sur- 
faces. In the red curve in Fig. [3] (b) for electron-doped 
graphene (n e = 0.04), the Kohn anomaly occurs both at 
qi between M and K and at q2 between K and T. The in- 
set in Fig.[3](b) shows the correspoding vectors qi and q2 
connecting the two electronic states at Fermi surfaces of 
two valleys in doped graphene. As the equibiaxial strain 
increases, the phonon at qi softens to zero frequency at 
the strain of e = 0.167 [blue line in Fig. |3] (b)]. This 
indicates that the doped graphene becomes structurally 
unstable through the less symmetric phonon mode at qi 
than A[ mode. The strain value of 0.167 for n e = 0.04 
doping differs from the corresponding critical breaking 
strain (0.179) obtained from the direct total energy cal- 
culation and stress-strain relationship [Fig. H](d)]. 

By taking stress and strain values when the softened 
phonon mode touches zero frequency, we can obtain the 
ideal strengths and the corresponding critical breaking 
strains as function of electron doping. As shown in 
Fig. [3] (c), the overall behaviors of ^Cidcai and <5e C ritical 
are similar to those obtained (Figure 2 (c) and (d)) from 
the stress-strain curves. However, the ideal strength in- 
creases only by ~ 1.7% and then decreases back to the 
original value as n e reaches 0.04 [Figure 3(c)]. The corre- 
sponding breaking strain shows the monotonic increase 
by ~ 15%. So, we can conclude that under equibiax- 
ial strain the electron and hole doping improve the ideal 
strength slightly or at least maintains its undoped value 
while the corresponding critical breaking strain increases 
dramatically by doping. These interesting behaviors of 
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(a) 



FIG. 4: (Color online) Stress-strain curves for graphene under 
(a) zigzag and (b) armchair strain with various doping levels. 
The doping levels are = 0.05, 0.02, undoped, n e = 0.02, 
0.05 from the top to the bottom for both plots, (c) The 
change of ideal strength and (d) critical breaking strain as a 
function of the doping level for the zigzag (open rectangles) 
and armchair (solid ones) strain respectively. 



stress and strain of doped graphene can be understood 
by noting that the ideal strength is mostly determined by 
a bonds between carbon atoms that does not be affected 
much by the doping but that the breaking point of strain 
is determined by 7r bands owing to the Kohn anomaly 
that changes greatly according to the doping. 

Next, we have analyized the ideal strength of doped 
graphene under uniaxial strain. Unlike equibiaxial strain 
case, the uniaxial strain breaks the crystal symmetry 
of graphene. Thus, under uniaxial strain, the pristine 
isotropic Dirac cone moves away from i^T-point and dis- 
torts into anisotropic form with direction dependent elec- 
tron group velocity 0, [To| ] . Here we only consider uniax- 
ial strains either along zigzag direction or armchair direc- 
tion The former and the latter are called as zigzag 
and armchair strain respectively. We calculate the stress- 
strain relationship to find out the doping dependent ideal 
strength under uniaxial strains. For these calculations, 
the primitive unit cell including two carbon atoms is used 
since uniaxially strained graphene breaks down elasti- 

Fig. I (a) and 



15] 



cally unlike equibiaxial strain cases 
(b) show the stress-strain curves for the zigzag and arm- 
chair strain respectively with different doping levels. It is 
found that the ideal strength under zigzag strain is sub- 
stantially more sensitive on the doping than that under 
the armchair one. Figs. 0] (c) and (d) show the change 
of ideal strength and critical strain as function of dopings 
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FIG. 5: (Color onlie) Electronic band structures of graphene 
under (a) the zigzag strain of e = 0.25 without doping, (b) 
e = 0.22 with electron doping n e = 0.05 and (c) e = 0.30 with 
hole doping nu = 0.05. The red and black lines are the a* 
and tt(tv*) bands respectively and the blue is the Fermi level. 
The inset in (a) shows the high symmetry points of BZ of 
uniaxially strained graphene. 



respectively. Contrary to the equibiaxial strain case, the 
response of graphene under uniaxial strain is asymmetric 
with respect to the electron and hole doping; graphene 
becomes stronger with the hole doping while it becomes 
weaker under the electron doping for both directions. For 
undoped graphene, the maximum stress of <7;deai = 40.2 
(36.6) occurs at e C riticai = 0.253 (0.198) for the zigzag 
(armchair strain) agreeing with a previous study [28j. 
As in the case for doped graphene under the equibiaxial 
strain, the doping changes the strain values a lot while 
the ideal strength changes a little upon doping [Figs. 4(c) 
and (d)]. Under zigzag strain, the critical breaking strain 
increases from 0.253 to 0.311 (23% increase) as the hole 
doping reaches = 0.05, and decreases to 0.22 (13% 
decrease) as the electron doping reaches n e = 0.05. On 
the other hand, under the armchair strain, it increases 
only by 3.5% and decreases by 5% in the same condition. 

Especially under the zigzag strain, the two Dirac cones 
at the K- and if'-point approach each other and eventu- 
ally merge after a critical strain 0, [13, I29I . This opens 
a small energy gap at the merged point |s], [l(| • At the 
same time, however, the energy of a* band is lowered 
toward the Fermi level [l(| signaling the breakdown of 
crystal structure of graphene. Our calculation shows that 
the stress maximum occurs as the unoccoupied a* band 
reaches the Fermi level. Figs. [5] (a), (b) and (c) show the 
electronic band structure for electron-doped, undoped, 
and hole-doped graphene under the zigzag strain respec- 
tively. For the electron doping with n e = 0.05 the a* 
band is closer to the Fermi level compared with one in 
undoped graphene under zigzag strain so that the band 
is filled up at a relatively lower strain of e = 0.22 [Fig. [5] 
(b)] . Contrary to that, for the hole doping with = 0.05 
the <7* band moves away from the Fermi level so that it 
requies relatively higher strain of e = 0.30 for the band to 
be filled up [Fig. [S](c)]. This explains the physical origin 
of electron-hole doping asymmetry of the ideal strength 
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of doped graphene under zigzag strain. 

Lastly, we discuss a possible solution to the aformen- 
tioned discre pan cy between the recent experiment [l3} 
and theories 1|1 1(| on the ideal strength of graphene. 
From our calculations, we note that the typical electron 
doped graphene can exhibit a theoretical ideal strength 
(> 40 N m _1 ) with a enhanced critical breaking cquib- 
iaxial strain, which is compatible with the experiment 
observation [L3[ . Though there is a still quantitative dif- 
ference in values of critical breaking strain between our 
theoretical estimation and experimental observation [l3| , 
the doping should play a major role to enhance the break- 
ing strain and other unknown experiment factors such as 
anisotropy in applied forces [1 31 ] or formation of disloca- 
tion 



14J may enhance the value further. 



In conclusion, we have studied the effect of doping on 
the mechanical properties of graphene. Electron- and 
hole-doped graphene is shown to improve or maintain 
its ideal strength under equibiaxial strain while the cor- 
responding breaking strain increases dramatically. So, 
it sets a new bound on the ideal strength and critical 
strain of isotropically expanded graphene in realistic sit- 
uations. Furthermore, we demonstrate that the uniaxi- 
ally strained graphene can decrease or increase the ideal 
strength and critical strain depending on the doping type. 
We believe that the present study not only highlights 
unique manifestation of electron and phonon physics of 
graphene but also can extend to explore interesting inter- 
play between electronic degree of freedom and mechanical 
properties of other two dimensional crystals [Io| • 
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